% Power
clc
clear all
%% subject 2 
cd ('E:')
load('Deltaband_d314posteo.mat')
X=Deltaband'; 
Fs=500; 
[Pxx,F] = periodogram(X,rectwin(length(X)),length(X),Fs);
save power_Deltaband_314posteo Pxx
p = bandpower(Pxx,F,'psd');
Postpower_Delta(1,2)=mean (p,'all');
Postpower_Delta(2,2)=min(Pxx,[],'all');
Postpower_Delta(3,2)=max(Pxx,[],'all');
Postpower_Delta(4,2)=median(Pxx,'all');
for b=1:19
p1=Pxx(:,b) ;     % building the vector
[pks1,locs1] = findpeaks(p1,F);
mpks1=find (pks1==max (pks1));
locs01=locs1(mpks1);
powerpeak(1,b) = locs01;
end
Postpower_Delta(5,2) = mean(powerpeak);

load Thetaband_d314posteo
X1=Thetaband';
Fs=500;
[Pxx1,F1] = periodogram(X1,rectwin(length(X1)),length(X1),Fs);
save power_Thetaband_d314posteo Pxx1
p01 = bandpower(Pxx1,F1,'psd');
Postpower_Theta(1,2)=mean (p01,'all');
Postpower_Theta(2,2)=min(Pxx1,[],'all');
Postpower_Theta(3,2)=max(Pxx1,[],'all');
Postpower_Theta(4,2)=median(Pxx1,'all');
for b=1:19
p12=Pxx1(:,b) ;     % building the vector
[pks1,locs1] = findpeaks(p12,F1);
mpks1=find (pks1==max (pks1));
locs01=locs1(mpks1);
powerpeak(1,b) = locs01;
end
Postpower_Theta(5,2) = mean(powerpeak);


load Alphaband_d314posteo
X3=Alphaband';
Fs=500;
[Pxx2,F3] = periodogram(X3,rectwin(length(X3)),length(X3),Fs);
save power_Alphaband_d314posteo Pxx2
p02 = bandpower(Pxx2,F3,'psd');
Postpower_Alpha(1,2)=mean (p02,'all');
Postpower_Alpha(2,2)=min(Pxx2,[],'all');
Postpower_Alpha(3,2)=max(Pxx2,[],'all');
Postpower_Alpha(4,2)=median(Pxx2,'all');
for b=1:19
p3=Pxx2(:,b) ;     % building the vector
[pks1,locs1] = findpeaks(p3,F3);
mpks1=find (pks1==max (pks1));
locs01=locs1(mpks1);
powerpeak(1,b) = locs01;
end
Postpower_Alpha(5,2) = mean(powerpeak);


load Betaband_d314posteo
X4=Betaband';
Fs=500;
[Pxx3,F4] = periodogram(X4,rectwin(length(X4)),length(X4),Fs);
save power_Betaband_314posteo Pxx3
 p03 = bandpower(Pxx3,F4,'psd');
Postpower_Beta(1,2)=mean (p03,'all');
Postpower_Beta(2,2)=min(Pxx3,[],'all');
Postpower_Beta(3,2)=max(Pxx3,[],'all');
Postpower_Beta(4,2)=median(Pxx3,'all');
for b=1:19
p4=Pxx3(:,b) ;     % building the vector
[pks1,locs1] = findpeaks(p4,F4);
mpks1=find (pks1==max (pks1));
locs01=locs1(mpks1);
powerpeak(1,b) = locs01;
end
Postpower_Beta(5,2) = mean(powerpeak);



load Gammaband_d314posteo
X5=Gammaband';
Fs=500;
[Pxx4,F5] = periodogram(X5,rectwin(length(X5)),length(X5),Fs);
save power_Gammaband_314posteo Pxx4
 p04 = bandpower(Pxx4,F5,'psd');
Postpower_Gamma(1,2)=mean (p04,'all');
Postpower_Gamma(2,2)=min(Pxx4,[],'all');
Postpower_Gamma(3,2)=max(Pxx4,[],'all');
Postpower_Gamma(4,2)=median(Pxx4,'all');
for b=1:19
p11=Pxx4(:,b) ;     % building the vector
[pks1,locs1] = findpeaks(p11,F5);
mpks1=find (pks1==max (pks1));
locs01=locs1(mpks1);
powerpeak(1,b) = locs01;
end
Postpower_Gamma(5,2) = mean(powerpeak);



load Deltaband_d314preeo
X6=Deltaband';
Fs=500;
[Pxx6,F6] = periodogram(X6,rectwin(length(X6)),length(X6),Fs);
save power_Deltaband_314preeo Pxx6
p06 = bandpower(Pxx6,F6,'psd');
Prepower_Delta(1,2)=mean (p06,'all');
Prepower_Delta(2,2)=min(Pxx6,[],'all');
Prepower_Delta(3,2)=max(Pxx6,[],'all');
Prepower_Delta(4,2)=median(Pxx6,'all');
for b=1:19
p6=Pxx6(:,b) ;     % building the vector
[pks1,locs1] = findpeaks(p6,F6);
mpks1=find (pks1==max (pks1));
locs01=locs1(mpks1);
powerpeak(1,b) = locs01;
end
Prepower_Delta(5,2) = mean(powerpeak);


load Thetaband_d314preeo
X7=Thetaband';
Fs=500;
[Pxx7,F7] = periodogram(X7,rectwin(length(X7)),length(X7),Fs);
save power_Thetaband_314preeo Pxx7
p7 = bandpower(Pxx7,F7,'psd');
Prepower_Theta(1,2)=mean (p7,'all');
Prepower_Theta(2,2)=min(Pxx7,[],'all');
Prepower_Theta(3,2)=max(Pxx7,[],'all');
Prepower_Theta(4,2)=median(Pxx7,'all');
for b=1:19
p07=Pxx7(:,b) ;     % building the vector
[pks1,locs1] = findpeaks(p07,F7);
mpks1=find (pks1==max (pks1));
locs01=locs1(mpks1);
powerpeak(1,b) = locs01;
end
Prepower_Theta(5,2) = mean(powerpeak);


load Alphaband_d314preeo
X8=Alphaband';
Fs=500;
[Pxx8,F8] = periodogram(X8,rectwin(length(X8)),length(X8),Fs);
save power_Alphaband_314preeo Pxx8
 p8 = bandpower(Pxx8,F8,'psd');
Prepower_Alpha(1,2)=mean (p8,'all');
Prepower_Alpha(2,2)=min(Pxx8,[],'all');
Prepower_Alpha(3,2)=max(Pxx8,[],'all');
Prepower_Alpha(4,2)=median(Pxx8,'all');
for b=1:19
p08=Pxx8(:,b) ;     % building the vector
[pks1,locs1] = findpeaks(p08,F8);
mpks1=find (pks1==max (pks1));
locs01=locs1(mpks1);
powerpeak(1,b) = locs01;
end
Prepower_Alpha(5,2) = mean(powerpeak);



load Betaband_d314preeo
X9=Betaband';
Fs=500;
[Pxx9,F9] = periodogram(X9,rectwin(length(X9)),length(X9),Fs);
save power_Betaband_314preeo Pxx9
 p9 = bandpower(Pxx9,F9,'psd');
Prepower_Beta(1,2)=mean (p9,'all');
Prepower_Beta(2,2)=min(Pxx9,[],'all');
Prepower_Beta(3,2)=max(Pxx9,[],'all');
Prepower_Beta(4,2)=median(Pxx9,'all');
for b=1:19
p09=Pxx9(:,b) ;     % building the vector
[pks1,locs1] = findpeaks(p09,F9);
mpks1=find (pks1==max (pks1));
locs01=locs1(mpks1);
powerpeak(1,b) = locs01;
end
Prepower_Beta(5,2) = mean(powerpeak);



load Gammaband_d314preeo
X10=Gammaband';
Fs=500;
[Pxx10,F10] = periodogram(X10,rectwin(length(X10)),length(X10),Fs);
save power_Gammaband_314preeo Pxx10
p10 = bandpower(Pxx10,F10,'psd');
Prepower_Gamma(1,2)=mean (p10,'all');
Prepower_Gamma(2,2)=min(Pxx10,[],'all');
Prepower_Gamma(3,2)=max(Pxx10,[],'all');
Prepower_Gamma(4,2)=median(Pxx10,'all');
for b=1:19
p010=Pxx10(:,b) ;     % building the vector
[pks1,locs1] = findpeaks(p010,F10);
mpks1=find (pks1==max (pks1));
locs01=locs1(mpks1);
powerpeak(1,b) = locs01;
end
Prepower_Gamma(5,2) = mean(powerpeak);
% 



Prepower_Beta(4,6)=median(Pxx9,'all');
for b=1:19
p09=Pxx9(:,b) ;     % building the vector
[pks1,locs1] = findpeaks(p09,F9);
mpks1=find (pks1==max (pks1));
locs01=locs1(mpks1);
powerpeak(1,b) = locs01;
end
Prepower_Beta(5,6)=mean(powerpeak);



load Gammaband_d326preeo
X10=Gammaband';
Fs=500;
[Pxx10,F10] = periodogram(X10,rectwin(length(X10)),length(X10),Fs);
save power_Gammaband_326preeo Pxx10
p10 = bandpower(Pxx10,F10,'psd');
Prepower_Gamma(1,6)=mean (p10,'all');
Prepower_Gamma(2,6)=min(Pxx10,[],'all');
Prepower_Gamma(3,6)=max(Pxx10,[],'all');
Prepower_Gamma(4,6)=median(Pxx10,'all');
for b=1:19
p010=Pxx10(:,b) ;     % building the vector
[pks1,locs1] = findpeaks(p010,F10);
mpks1=find (pks1==max (pks1));
locs01=locs1(mpks1);
powerpeak(1,b) = locs01;
end
Prepower_Gamma(5,6)=mean(powerpeak);



